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PREFACE 

This  report  examines  a technique  for  measuring  rocket  motor  thrust 
malalignment  by  observing  the  associated  motion  of  a tethered  (at  a 
single  point)  rocket.  The  technique  was  proposed  by  Mr.  Tommy  Howell 
of  the  Test  and  Evaluation  Directorate,  U.  S.  Army  Missile  Research, 
Development  and  Engineering  Laboratory  at  Redstone  Arsenal.  The  analysis 
reported  herein  has  been  supported  by  that  Directorate  through  an  Army 
Research  Office  grant. 

The  method  examined  here  was  proposed  so  as  to  utilize  rotational 
dynamics  as  a response  mechanism  instead  of  more  traditional  methods 
of  suspension  for  testing  which  usually  yield  a combination  of  motions. 
Indeed,  most  previous  systems  required  a six-degree-of  freedom  analysis 
to  relate  observed  motion  to  thrust  malalignment. 

We  gratefully  acknowledge  the  assistance  of  Messrs.  Tommy  Howell, 

Ted  Haas  and  Lei  and  Lee  of  the  Redstone  Arsenal  who  have  provided 
technical  data  for  this  study. 

We  also  acknowledge  the  editorial  skill  and  extreme  good  humor  of 
Mrs.  Dianne  Fretwell  and  Mr.  Walter  Clement  of  the  Auburn  University 
Mechanical  Engineering  Department  who  were  indispensable  in  preparing 
the  manuscript. 

Project  Leader 


Joe  W.  Reece 
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ABSTRACT 

The  lateral  perturbine  forces  due  to  thrust  malalignment  can 
significantly  affect  the  flight  characteristics  of  a missile.  Past 
attempts  to  predict  thrust  malalignment  have  not  proven  successful. 

This  report  proposes  a new  method  for  making  valid  measurements  of 
the  malalignment  phenomena  through  the  use  of  a single  point  tethered 
missile  technique.  The  equations  of  motion  for  the  missile  in  the  test 
device  are  developed,  coded  into  a computer  and  tested  against  well 
known  exact  analytical  solutions  for  motion  about  a point.  The  mass, 
inertia  and  thrust  equations  for  a typical  test  missile  are  then  coded 
into  the  computer  and  the  resultant  motion  paths  are  simulated  for 
various  initial  spin  rates  and  degrees  of  thrust  malalignment.  Analysis 
of  the  ccllected  data  indicates  that  the  tethered  technique  is  a viable 
solution  for  predicting  thrust  malalignment. 
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I.  INTRODUCTION 

A major  problem  in  missile  dynamics  is  unpredictable  flight  path 
motion  due  to  lateral  perturbing  forces.  Thrust  vector  malalignment, 
which  occurs  when  a missile's  thrust  vector  is  not  coincidental  with 
its'  longitudinal  spin  axis,  is  often  a cause  of  these  forces.  In  order 
to  reduce  the  effects  of  thrust  malalignment,  spin  stabilization  is  often 
employed.  The  missile  is  given  a spin  rate  about  its'  longitudinal  axis 
which  rotates  the  various  perturbing  forces  and  makes  them  less  signi- 
ficant. 

Since  the  effects  of  thrust  malalignment  are  a serious  problem  in 
achieving  a specified  missile  flight  path,  it  is  desirable  to  be  able 
to  measure  the  degree  of  malalignment  in  a particular  missile  system. 

The  knowledge  thus  gained  can  be  used  to  correct  design  characteristics 
and  produce  better  downrange  flight  motions. 

Various  attempts  to  measure  thrust  malalignment  have  been  made. 

The  techniques  which  have  been  used  in  the  past  have  not  produced  use- 
ful results.  A new  proposed  method  for  measuring  thrust  malalignment 
is  indicated  in  Figure  1.  A missile  for  testing  is  suspended  in  a ball 
and  socket  joint  which  allows  three  degrees  of  freedom.  The  spin  and 
main  motors  of  the  missile  are  started  sequentially.  The  missile  may 
be  constrained  initially  as  desired  to  give  uniform  starting  conditions. 
Depending  on  the  amount  of  malalignment  and  initial  spin,  the  missile 
will  undergo  some  particular  subsequent  motion.  A circular  constraint 
ring  is  provided  to  prevent  excessive  lateral  excursions. 
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The  objective  of  this  report  is  to  predict  the  motion  of  a test 
missile,  in  the  constraining  fixture,  with  known  inertia,  mass,  and 
thrust  characteristics  for  various  degrees  of  thrust  malalignment  and 
initial  spin  rates.  The  parameters  describing  the  motion  of  the  missile 
are  displayed  at  the  end  and  may  serve  as  a set  of  prediction  curves  for 
actual  testing.  If  one  measures  the  motion  of  a real  missile  in  the 
test  fixture,  the  curves  may  be  used  to  relate  the  degree  of  thrust 
malalignment.  The  results  given  here  can  also  be  used  as  a basis  for 
further  consideration  for  this  method  of  testing. 

The  results  of  this  analysis  show  that  this  testing  method  may  be 
a viable  way  to  measure  thrust  malalignment.  Excursions  predicted  for 
a typical  rocket  are  within  the  measurable  range  and  do  not  take  the 
rocket  into  unreasonable  motion  patterns.  Also,  some  surprising  results 
occur  which  shed  some  light  on  the  necessary  spin  rate  required  for  the 
stabilization  of  rockets. 

The  analysis  involves  five  steps: 

1.  Development  of  the  equations  of  motion  for  the 
constrained  missile. 

2.  Development  of  a technique  for  solving  the  equations 
of  motion. 

3.  Testing  of  the  solution  technique  against  known 
exact  solutions  (Appendix  II). 

4.  Establishment  of  the  basic  characteristics  for  a 


typical  test  missile.  These  characteristics  will 
include  mass,  inertia,  and  thrust  equations  as  a 
function  of  time. 


4 


5.  Coding  into  a computer  the  equations  in  (4)  for  the  test 
missile  to  generate  a set  of  curves  for  predicting  its' 
motion  in  the  test  fixture. 


II.  EQUATIONS  OF  MOTION 


The  motion  of  a missile  about  a fixed  point  can  be  described  by 
three  Euler  angles,  e,  (t>,  [1].  These  three  angles  are  used  to  give 
the  spatial  interrelationship  between  three  coordinate  systems,  one  of 
which  is  attached  to  the  missle.  Figure  2 shows  the  three  Euler  angles 
and  their  relationship  to  the  three  coordinate  systems. 

The  primary  inertia  reference  frame  will  be  the  XYZ  coordinate 
system.  It  is  the  reference  frame  from  which  an  observer  would  record 
the  motion  of  the  missile  relative  to  the  earth.  A second  reference 
frame,  the  x'y'z'  set,  is  rigidly  attached  to  the  missile.  The  forces 
and  associated  moments  along  with  the  mass  and  inertia  properties  of  the 
missile  are  described  by  this  reference  frame.  Finally,  an  intermediate 
coordinate  system,  the  xyz  set,  is  used  to  relate  the  location  of  the 
missile  longitudinal  spin  axis  relative  to  the  XYZ  inertia  reference 
frame.  Each  of  the  reference  frames  have  a common  origin  at  0. 

The  xyz  reference  frame  is  related  to  the  XYZ  inertia  frame  by  the 
two  Euler  angles  ii  and  9.  In  general,  coordinate  axis  x of  the  xyz  set 
moves  in  the  XY  plane.  Its  angular  relationship  to  the  X axis  is  given 
by  the  precession  angle  ip.  Also,  the  y and  z axis  of  the  xyz  set  are 
allowed  to  rotate  about  the  x axis.  The  angular  relationship  between 
the  Z and  z axes  is  given  by  the  nutation  angle  e. 

Coordinate  system  x'y'z'  may  be  described  by  a single  Euler  angle 
(j)  relative  to  the  xyz  frame,  since  the  z and  z'  axes  are  always  coin- 
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cidental.  The  angular  relationship  between  the  x and  x'  axes  is  given 
by  the  spin  angle  <t). 

Associated  with  each  of  the  Euler  angles  i;,  e,  (t>  are  the  angular 
velocities  ip,  e,  41  and  the  angular  accelerations  p,  e,  <p.  The  nose  of 
the  missile  is  considered  fixed  at  the  origin  0 and  its'  longitudinal 
spin  axis  lies  along  the  -2'  coordinate  axis  as  in  Figure  2. 

For  the  coordinate  systems  described,  the  degree  of  missile  thrust 
malalignment  is  given  by  the  angle  a [rad.],  which  is  the  angular  measure 
between  the  missile's  thrust  vector  and  its'  longitudinal  spin  axis. 

The  total  force  acting  on  the  missile  at  any  time  is  given  by  eF. 

Moments  due  to  the  various  forces  are  taken  about  the  origin  0 and  the 
total  moment  is  given  by  eM^.  The  direction  of  the  force  due  to  gravity 
is  taken  along  the  -Z  axis  for  all  time. 

In  general,  the  angular  momentum  of  the  missile  is  given  by  the 
realtion  [2]: 

ff  = / ^ Vdm  ( 1 ) 

b 

where, 

^ = position  vector  to  a differential  mass  dm 

V H velocity  vector  of  the  mass  dm. 

Since  the  missile  maintains  symnetry  with  respect  to  its'  mass  distribu- 
tion along  the  z'  axis,  a principle  axes  problem  exist  and  therefore, 
the  products  of  inertia  vanish.  The  angular  velocity  vectors  n and  u 
may  be  used  to  describe  the  angular  motions  of  the  xyz  and  x'y'z'  sets, 
respectively,  relative  to  the  XYZ  coordinate  set,  as  shown  in  Figure  3. 
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The  angular  velocity  of  the  xyz  set  is  given  by  [1]: 


n = 0e+i(/sin0e  + \p  cos  0 e 
X y z 


and  the  angular  velocity  of  the  x'y'z'  set  is  given  by: 


u)  = 0 e + 4)  sin  e e + (4)  + 4*  cos  e)  e 
X y z 


Since  the  velocity  V in  equation  (1)  is  specified  by: 


V = w X 


and  there  are  no  products  of  inertia,  then. 


x'x'  X X y y y y z'z  z z 


The  inertia  components  are  therefore  given  as. 


Ix-x-  = J (Py-'  + PZ-')  d"! 


Iy>y'  = J (p^'^  Px'^)  dm 


III  = / (p„i^  + Pw'^)  dm 


(2) 


(3) 


(4) 

(5) 

(6) 


Since,  syimetry  along  the  z' 


exists,  then: 

'zz  = 

'z-z- 

= I 

^x  - 

I , 

XX 

= I 

I 1 1 
y y 

= I 

axis  with  respect  to  mass  distribution 
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and, 

= re'  (7) 

Hy  = Tij/  sin  9 (8) 

= I ((}i  + i//  cos  e)  (9) 

The  derivative  of  the  momentum  components  may  be  expressed  as: 

= i'e  + re  (10) 

Hy  = I'li/  sin  9 + I'i|;  sin  e + I'lf/  e cos  e (11) 

= I (41  + cos  e)  + I (<))  + 0 cos  e - \p  e sin  e)  (12) 

The  general  equation  of  motion  for  the  body  is  given  as  [2]: 


m=^^t  +nxH 
^ xyz 


(13) 


The  scalar  components  of  the  moment  vector  equations  given  by  (13)  are: 


= I'e  + I'e  + Lj>ip  sin  e + (I-I')  sin  e cos  e (14) 

My  = I'ii(  sin  9 + I'\p  sin  e + Zl'ij^e  cos  9 - Ie(4i  + cos  e)  (15) 


M_  = I 4>  + U + U cos  9 - hpe  sin  e + I ip  cos  e 


(16) 


where  M^,  M^,  M^  are  the  components  of  zM^.  The  angular  accelerations 
are  then: 


M^  - I'e  - sin  e - (l-I')  4;^sin  9 cos  e 


(17) 


M - I'll)  sin  9 - 21' 1// 9 cos  e + le  (41  + 0 cos  e) 

I*  = 


r sin  e 


(18) 


n 


M - 1 4,  - I if)  cos  0 - I ij)  cos  0 + I if)  0 sin  0 

4=-^ (19) 

In  general,  the  moments  in  the  x'y'z'  body  fixed  coordinate  system  are 
known.  The  x'y'z'  moment  components  can  be  transformed  into  moments  in 
the  xyz  set  by  the  matrix  operation: 


M 

X 

M 

_ 

y 

M, 

z 

cos 

- sin 

4’ 

I 

0 

sin 

cos 

4> 

0 

• 

1 

0 

0 

1 

_«z'_ 

which  results  in  the  relations: 


M = M , cos  <j)  - M , sin  41  + wR  sin  0 

A /\  y 

My  = Sin  4>  + M^,  cos  <p 


If  the  following  variables  are  defined: 


V = if)  cos  0 
A = if)  sin  0 
e = cos  0 
y = sin  0 
e = cos  4> 
p = sin  4> 


(21) 

(22) 

(23) 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


Equations  (21)  through  (30)  and  equations  (17),  (18),  and  (19)  may 
be  simplified  to: 
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(31) 

(32) 

(33) 


A suitable  integration  technique  is  required  to  solve  these  nonlinear 
differential  equations  (31)  through  (33). 

To  predict  the  location  of  a fixed  point  along  the  spin  axis  of 
the  missile,  the  relations: 


III.  INTEGRATION  TECHNIQUE  FOR  SOLVING  THE  EQUATIONS  OF  MOTION 


The  solution  algorithm  for  equations  (31),  (32),  and  (33)  involves 
the  use  of  a numerical  forward  marching  technique.  In  general,  the 
angular  positions,  velocities  and  accelerations  of  the  various  coordinate 
systems  are  time  dependent.  Therefore,  their  components  may  be  expressed 


'I'  = f^  (t)  ; 'll  = f2  (t)  ; = ■•'3 

9 = g,  (t)  ; 9 = gp  (t)  ; 9 = g_  (t) 


h^  (t)  ; ({)  = h^  (t)  ; 


= h3  (t) 


Consider  the  precession  acceleration  and  velocity  which  may  be  expressed 


At-*-o  At 


Aij<  (t)  ^ Lim  ^^2^  ~ ^^1^ 


At  -»■  0 


■ _ Lim  A0  (t)  _ Lim 


Ip  (^2^  ~ 'P  (T-j) 


At->-0  At 


At  0 


where  At  = t2  - . 

Equations  (36)  and  (37)  may  be  rewritten  as; 


Lim  i|;  At  = 0 (tp)  - \p  (t-,) 
At  -*■  0 


Lim  tj;  At  = 4-  (tp)  = Ip  (t,) 
At  > 0 ^ 


where  4;  and  4/  may  be  evaluated  at  t = t-j 
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As  the  time  interval  At  decreases,  the  approximations  given  by: 

1//^  M = i!  [t^)  - Ip  (t^)  (40) 

At  = i|)  (t^)  - Ip  (t^)  (41) 

become  more  exact. 

Equations  such  as  (40)  and  (41)  may  be  coded  for  computed  solutions 
to  equations  (31),  (32),  and  (33).  The  basic  relations  for  the  computer 


are: 

\p  (12)  = At  + (t-| ) (42) 

Ip  (t2)  =1^1  At  + i|/  (t-| ) (43) 

0 (t2)  =0^  At  + 0 (t-| ) (44) 

0 (t2)  =6-1  At  + 0 (t^ ) (45) 

(J)  (t2)  = (61  At  + (Ji  (t-j)  (46) 

(t2)  = At  + (^  (t^)  (47) 


Appendix  I gives  the  computer  code  used  to  solve  equations  (31)-(33) 
by  equations  (42)-(47).  The  basic  solution  may  be  outlined  as 
follows: 

(1)  Establish  initial  missile  parameters  such  as  orientation, 
mass  and  inertia  properties.  Also  establish  t = 0 as  the  solution 
starting  point  in  time. 

(2)  Given  the  existing  set  of  spatial  constraints,  calculate 


the  moments  acting  upon  the  missile  about  the  origin  0. 
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(3)  Solve  for  the  new  angular  acceleration  components  of  the 
reference  frames  using  equations  (31),  (32),  and  (33). 

(4)  Solve  for  the  new  angular  position  components  of  the  reference 
frams  using  equations  (43),  (45),  and  (47). 

(5)  Solve  for  the  new  angular  velocity  components  of  the  reference 
frames  using  equations  (42),  (44),  and  (46). 

(6)  Increment  the  time  variable  of  the  solution  by  aT  and  check 
to  see  if  the  elapsed  time  is  greater  than  the  time  requested  for  the 
solution.  If  the  elapsed  time  is  less  than  the  time  requested,  go  to 
step  (2).  Also,  print  out  any  results  obtained.  If  the  elapsed  time 
has  exceeded  the  time  requested  for  th'^  solution,  end  the  program. 
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IV.  THRUST,  MASS  AND  INERTIA  PARAMETERS 
FOR  A TYPICAL  TEST  MISSILE 

For  investigating  the  motion  of  a typical  missile  system,  a set 
of  test  parameters  for  thrust,  mass  and  inertia  properties  have  been 
developed.  Figure  5 shows  the  thrust  versus  time  relation  for  a missile 
under  investigation.  The  thrust  simulation  consists  of  four  linear 
segments  which  comprise  four  special  time  phases  of  the  burning  process. 


These  relations  may  be  expressed  as  follows: 

Thrust  = SOOOOOt  (lb)  0 ^ t £ .05  sec  (48) 
Thrust  = 18333. 33t  + 14000  (lb)  .05  £ t £ .65  sec  (49) 
Thrust  = 26000  (lb)  .65  < t < .90  sec  (50) 
Thrust  = -520000t  + 494000  (lb)  .90  £ t < .95  sec  (51) 
Thrust  = 0 (lb)  t > .95  sec  (52) 


The  time  dependent  thrust  relations  can  be  used  to  generate  a set 
of  moments  about  the  origin  0 in  the  x'y'z'  coordinate  set.  If  a is 
defined  as  the  thrust  malalignment  angle,  measured  in  radians,  then 
the  moment  relations  acting  upon  the  missile  may  be  given  as: 


(53) 

= 0 (in-lb) 

(54) 

= 0 (in-lb) 

(55) 
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Thrust  Versus  Time  Relation  for  a Typical  Missile 
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For  the  missile  under  investigation,  = -^35  in.,  is  defined 

as  the  location  on  the  missile  where  the  thrust  vector  is  considered  to 
act  as  a point  force.  Therefore,  a set  of  moment  versus  time  relations 
can  be  developed  by  substituting  equations  (48)  through  (52)  into 
equation  (53): 

= 40500000t  sin(ll .538ta)  (in-lb)  0 < t < .05  sec  (56) 

M^,  = (2475000t  + 1890000)sin((.7051t  + .5385)a)  (in-lb) 

.05  ^ t £ .65  sec  (57) 

M^,  = 3510000  sin  (a)  (in-lb)  .65  < t < .90  sec  (58) 

M^,  = (-70200000t  + 66690000)  sin  ((-20t  + 19)a)  (in-lb) 

.90  £ t £ .95  sec  (59) 

= 0 (in-lb)  t ^ .95  sec  (60) 

The  weight  of  the  missile  varies  because  of  fuel  burn-up  as, 

w = -87.589t  + 250.19  (lbs)  (61) 

and  the  center  of  gravity  location  in  the  x'y'z'  set  can  be  expressed 
as: 

= 14.3011t  - 61.518  (in)  (62) 

Therefore,  the  moment  due  to  gravity,  M = Rg2”>  varies  according  to 
the  equations, 

M = -1252.619t2  + 8966. 3t  - 15391.2  (in-lb) 

0 £ t ^ .95  sec  (63) 


M = -8003.7  (in-lb) 


t > .95  sec 


(64) 
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The  roll  inertia  of  the  missile  is 

Voll  ” -14-614t  + 43.155  (slug-in^)  0 ^ t ^ .95  sec  (65)  ' 

= 29.271  (slug-in2)  t > .95  sec  (66) 

While  the  transverse  inertia  varies  according  to  the  equations, 

^trans  ~ -22428. 59t  + 39196.2  (slug-in^)  0 ^ t .95  sec  (67) 

i 

Itrans.  ^ '•7889.75  (slug-in^)  t > .95  sec  (68)  j 

Equations  (56)  through  (60),  (63)  through  (68)  were  coded  into  the 
computer  to  produce  the  results  shown  in  section  V. 

I 

I 

I 

i. 

! 

i 

I 

1 


1 
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V.  RESULTS  OF  TEST  MISSILE  MOTION  STUDIES 


Figures  6 through  23  are  the  result  of  coding  the  data  in  section 
IV  into  the  computer.  Four  initial  spin  rates  of  0,  5,  10,  and  15  rps 
(rev/sec)  were  utilized  in  making  the  plots.  For  each  initial  spin 
rate,  five  values  of  maximum  thrust  malalignment,  .001,  .002,  .003,  .004, 
and  .005  radians  were  used. 

In  the  case  of  ij)  = 0 rad/sec  initial  spin  rate,  the  nutation 
components  consisting  of  e,  e,  e were  plotted  versus  time.  For  the 
initial  spin  rates  of  a = 5,  10,  and  15  rps,  the  position  of  the  missile 
nozzle  located  at  z'  = -135  in.  was  plotted  in  XY  coordinates  with  the 
functional  parameter  of  time  indicating  the  starting  and  ending  points 
of  the  solution.  An  observer  located  at  the  origin-0,-see  Figure  2, 
and  recording  each  event  would  notice  the  given  nozzle  patterns  in  the 
z-direction. 

Test  (a)  Thrust  Malalignment  of  a Test  Missile  with  no  Initial 
Spin,  !|i  = 0 rev/sec. 

Figures  6,  7,  and  3 give  the  nutation  angle,  velocity,  and 
acceleration  of  the  test  missile  without  an  initial  spin  rate  (41  = 0 rps). 
The  various  stages  of  the  thrust  process  are  easily  observed  in  the 
nutration  acceleration,  (see  Figure  8).  As  the  rate  at  which  thrust 
is  generated  changes  abruptly,  discontinuities  in  the  nutation  accelera- 
tion curves  occur.  Specifically,  at  times  t = .05  sec,  t = .65  sec, 
t = .90  sec,  and  t = .95  sec,  abrupt  changes  in  the  thrust  generation 
rate  result  in  discontinuous  changes  in  the  nutation  acceleration 
curves.  Observation  of  the  nutation  angle  versus  time  curve  (see 
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Fiqure  6.  Nutation  Angle  as  a Function  of  Time  (41  = 0 rad/sec) 


Figure  7.  Nutation  Velocity  as  a Function  of  Time  (f  = 0 rad/sec) 


Figure  9.  Nozzle 


COORD 
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Figure  6),  indicates  that  large  deflections  of  the  missile  will  occur. 

Typically,  deflections  on  the  order  of  .1  radian  will  be  observed.  For 
anticipated  malalignment  values  of  a = .005  radian,  a constraint  ring 
may  be  required  to  prevent  excessive  deflections  of  the  missile  in  the 
test  fixture. 

Test  (b)  Thrust  Malalignment  of  a Test  Missile  with  5 rps  Initial 
Spin  Rate  ((f,  = 5 rev/sec). 

Figures  9 through  13  indicate  the  paths  of  motion  of  the  nozzle  of 
the  constrained  missile.  The  nozzle  is  considered  to  be  located  at  the 
point  -135  in.  along  the  z'  axis  from  the  origin  0.  After  the  main 
motor  is  fired,  the  thrust  malalignment  of  the  missile  results  in  a 
lateral  perturbing  force  which  rotates  about  the  spin  axis.  This  force 
generates  a moment  which  when  coupled  with  a rather  large  moment  due  to 
gravity,  (equation  63),  produces  a spiralling  path  of  motion. 

This  spiralling  path  of  motion  tends  to  deflect  outward  due  to  the 

j 

combined  moments  of  gravity  and  thrust.  The  effects  of  spin  stablilza- 

tion  are  well  demonstrated  in  the  plots  since  the  nutation  of  the  missile 

is  significantly  reduced.  Also,  four  spirals  are  observed  in  each  of  the  i 

I 

plots  which  represent  the  effect  of  rotating  the  lateral  forces  five 
times  before  motor  burnout. 

Note  that  the  cumulative  effect  is  to  move  the  nozzle  in  one 
direction,  i.e.,  the  motion  is  not  centered  about  the  initial  position 
of  the  missile.  This  leads  to  the  conjecture  that  the  spin  rate  may 
be  far  too  low  to  negate  the  effects  of  thrust  malalignment  for  a 
rocket  motor  which  burns  out  so  quickly.  This  observation  is  applicable 
to  practically  all  the  test  investigated  in  this  work.  Of  course,  the 
extremely  short  burn  time  of  .95  sec.  should  also  be  recalled. 

__ J 


MOTOR  FIRED  (t=0  sec) 


.002  rad  THRUST  MALALIGNMENT 
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Figure  15.  Nozzle  Position  as  a Function  of  Time  (10-.002) 


Figure  16.  Nozzle  Position  as  a Function  of  Time  (10-.003) 


Fiqure  17.  Nozzle  Position  as  a Function  of  Time  (10-.004) 


MOTOR  FIRED  (t=0  sec) 


ROINfll 
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Test  (c)  Thrust  MalaUgnnient  of  a Test  Missile  with  10  rps  Initial 
Spin  Rate  (<t>  • 10  rev/sec). 

Figures  14  through  18  simulate  the  motion  path  of  the  missile  nozzle 
located  at  -135  In.  along  the  z'  axis  from  the  origin  0.  As  stated  before, 
the  spiralling  motion  Is  due  to  a combined  lateral  thrust  moment  and  a 
gravitational  moment.  Several  Interesting  points  may  be  noted  about  the 
figures.  The  most  obvious  effect  observed  is  the  significant  reduction 
In  nutation  as  a result  of  a higher  spin  rate.  Also,  the  radial  path  of 
deflection  has  rotated  counterclockwise  from  the  path  of  deflection  for 
the  5 rps  case.  A total  of  eight  spirals  are  observed  In  the  figures. 

Each  of  these  effects  may  be  attributed  to  t(ie  incre<jscd  spin  ra*e  of 
tbe  constrained  missile 

lest  (d)  ThtiiSt  Malal  ignmer.f-  of  a lest  Missile  with  15  rps  Ir-itial 
Spin  Rate  (i>  = 15  rev/sec) . 

Figures  ’9  through  23  show  ttie  path  of  motion  for  the  nr/zle 
located  as  in  ‘est  ( asos  (b)  ard  'c).  Ihe  increased  '.pin  ate  le  ults 
in  further  counterclockwise  rotation  of  t!ie  tadial  cat')  of  deflection 
for  the  spiralling  motion.  Also,  the  numl  er  of  complettd  spi-als  has 
increased  to  thirteen  in  each  of  the  plots.  Iwo  observations  mav  be 
made  which  are  the  result  of  a higier  spin  rate.  First,  there  is  a 
significant  further  reduction  in  nutation  below  that  found  in  cases 
(b)  and  (c)  Second,  the  spiralling  path  of  motion  tends  to  cross 
itself  much  more  since  the  effects  cf  the  lateral  perturbi'o  forces 
have  been  decreased. 

lest  (e)  Results  from  a Square  Wave  Ihrust  Impulse. 

The  purpose  of  this  test  was  to  study  the  motion  of  a missile  with 
identical  inertia  and  mass  characteristics  as  ir  cases  (a)  through  (d). 


.001  rad  THRUST  MALALIGNMENT 


Figure  19.  Nozzle  Position  as  a Function  of  Time  (15-.001) 


.002  rad  THRUST  MALALIGNMENT 


Figure  20.  Nozzle 
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except  that  the  thrust  relation  is  a square  wave.  Consider  a missile 
with  a thrust  output. 

Thrust  = 26000  lb  ^ ^ -50  sec 

Thrust  = 0 lb  t > .90  sec 

The  moments  generated  in  the  x'y'z'  set  are  given  as: 


M . 

X 

= 3510000  sin  (a)  (in-lb) 

0 < t < .90  sec 

(69) 

M , 

X 

= 0 (in-lb) 

t > .90  sec 

(70) 

M , 

y 

= 0 

= 0 

Equations  (63)  through  (70)  were  coded  into  the  computer  along  with  an 
initial  spin  rate  of  $ = 5 rps  for  a = .001 -.005  radians  thrust  malalign- 
ment. The  results  are  shown  in  Figures  24  through  28  and  should  be 
interpreted  as  in  Figures  9 through  13.  A comparison  of  these  figures 
indicates  the  importance  of  knowing  the  thrust  versus  time  relationship 
of  the  missile.  If  an  actual  test  is  to  be  performed,  it  may  be  necessary 
to  monitor  the  thrust  output  as  a function  of  time  to  get  accurate 
results.  The  comparison  of  the  two  sets  of  figures  show  radical  dif- 
ferences in  the  nozzle  oaths  of  motion.  This  discrepancy  is  the  result 
of  the  differences  in  the  thrust  versus  time  relations  for  the  two  cases. 


Figure  ?A . Nozzl  Position  as  a Function  of  Time  (S5-.001) 


Figure  25.  Nozzle  Position  as  a Function  of  Time  (S5-.002) 


Figure  27.  , Nozzle  Position  as  a Function  of  Time  (S5-.004) 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Results  from  the  simulation  studies  for  a test  missile  indicate 
that  maximum  nutation  occurs  when  the  initial  spin  rate  is  zero.  If 
the  maximum  thrust  malalignment  is  expected  to  exceed  5 mrad,  then  a 
constraint  ring,  (see  Figure  1),  should  be  considered  for  the  test 
fixture.  Results  from  the  test  cases  also  indicate  that  a knowledge  of 
the  thrust-time  relationship  is  essential  for  making  accurate  predictions 
of  the  results  of  thrust  malalignment.  Also,  the  inertia  and  mass 
characteristics  for  the  missile  during  the  motor  burning  phase  are 
required. 

The  most  accurate  way  of  predicting  thrust  malalignment  appears  to 
be  the  most  simple.  If  the  missile  is  not  given  an  initial  spin,  maxi- 
mum nutation  will  occur.  Therefore,  a very  precise  measurement  of  mal- 
alignment may  be  made.  This  study  neglected  all  aerodynamic  forces  and 
frictional  moments  at  the  pivot  joint  of  the  test  fixture.  If  the 
actual  constraining  device  is  built,  consideration  of  these  forces 
and  moments  may  be  included  easily  in  the  moment  equations  for  the 
missile. 

It  should  be  recalled  that  the  objective  of  this  investigation  was 
to  determine  the  range  of  motion  resulting  from  anticipated  thrust 
malalignment  anomalies  affecting  a typical  rocket  suspended  in  this 
unique  ball  and  socket  fashion.  The  results  of  the  analysis  show  that: 

• The  motion  ranges  and  motion  rates  are  in  the  measurable 
regime. 
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• This  scheme  may  be  a viable  experimental  method  for 
gathering  performance  data  on  real  firings. 

Further  work  needs  to  be  done  to  examine  questions  of  the  effects 

t Combined  Unbalance  and  Thrust  Malalignment 
t Vibration  and  Buckling  Phenomena 

• Hardware  Imperfections  in  the  Mounting  and  Suspension  Rig 
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Thrust  Malalignment  Computer  Code 
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APPENDIX  I 


Thrust  Malalignment  Computer  Code 

The  computer  code  for  predicting  the  resultant  motion  of  a missile 
for  various  degrees  of  thrust  malalignment  and  initial  conditions 
utilizes  the  following  list  of  vanaoles.  The  computer  variable  is 
listed  first,  followed  by  the  mathrndtical  equivalent. 

Ml  = M , 

X 

M2  = My,  moments  in  the  x'y'z'  coordinate  system 

M3  = M^. 

M4  = M I component  in  the  xyz  coordinate  system 

II  = I I 

^x'z'  component  and  derivative  (roll  inertia) 

^x'x'  ~ ^v'v'  component  and  derivative 

^ ^ (Transverse  Inertia) 


A1  = «|/ 

A2  = 

A3  = '■!; 

j I number  of  aT  intervals  between  successive  printings 

of  output 

I j total  number  of  output  printings  ] 
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Euler  angle  precession  components 


13  = r 

14  = i' 


I 
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A4  = 0 
A5  = e 
A6  = 0 


Euler  angle  nutation  components 


A7  = 

A8  = (^ 
A9  = ^ 


Euler  angle  spin  components 


B1  = 0 
B2  = li 
B3  - V 
B4  = X 
B5  = e 
B6  = p 


basic  control  variables 

(see  equations  (24)  through  (29)) 


X = X 

Y = y 


location  of  the  missile  nozzle  in  the  XYZ 
coordinate  system 


R = R 


-z  location  of  the  missile  nozzle 


A = a 


thrust  malalignment  (radians) 


T = at 


integration  time  inverval 


T1 


total  elapsed  time 
Note:  T1  = I*J*T 
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I R ( 11 . L T . . < , ) ft  1 = 1 j 1 : j . « s I \ ( A ) 

I F ( r I .LT  . .S  J ) ft  1=  (2^7-3CJ2.*r  1 ♦ISS:::C>.  . I’J'S  I*7(  ( . 7s  51  *n  ♦.  -3  'P‘j  )«A  ) 

IF  ( T 1 .LT  . . 3 ) ''I  = A . S5E-*-7*T  1»5  livil  1 I .5  39>!=r  1 #A  ) 

c l'NL  CF  SRAllAL  CCNSTRAINTS 

-•l  = Cf  S(ft'0 
b'  2 = S I ft  ( \ A ) 

I-  .1  = A 2 ft  P 1 
t-<  = A2*fV 
'•  -3  = c r.  s ( A 7 ) 
r 6 = S I \ ( A 7 ) 

M'  ( f ft  . ' S . , . ) . r T r 1 

AT=(''lftl'iSft''2ft.  5-UftEM+'  I lft(AP  + HJ)-2.ftI  UP3»ftfl5)/( 

I 7.  - ( ft  1 ftps  - ft?ft  r t ift  AftP2  - I A>i>A5  + { 1 Jfth  T - I 1ft  ( AH4  P T M ftPA  ) / ( T 

A',  = ( ft  1-  1 ^ » I AP+P  3 ) - lift  { A 3*P  l-AbftHvi  ) ) / I 1 


i 
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<51  = i?*  r ♦ A I 
A2  = A1<‘T  + A2 
+ A^ 

A5  = A6^T  + A’j 
A7-A{^«I  + A7 
Ae  = <,'5<>r  + AB 

2 T 1 = T 1 ♦ r 

I Tf  ( U , 'i  ) T 1 , AA  , A5  , A6 

3 1 C K iv  A 1 ( ^ F,  1 5 . A ) 

A C(.  M IM'c 

STCP 

FM 

//CL.FI  ICFCCl  n:  LS;\  = ^ECAfi. BINGO, Cl  SP  = (N£A,C4rLC)  , 

//  L ApcL  = ( 1 , SL  ) ,LN  I T = T4P£  ,CCP=  ( RtCF>=(-  RS  , LPI-a  = 6:^  ,H  Cl  / £=  1 AC2.,  ) 
//Gc.svsiN  rc  <■ 

//STEP2  FXtfC  T".  PEFU.E 
//CYSIN  CC  » 

FILE  r SN  am-  =Rr  C^iP  . T’  IFGC. 


t?jCf  CPITIYE=  L.12SFC 


//<l.;.K'T  JLi'  ( ‘'r  L'.H  , FC'"«'S  = 6^Sri  ) , • Chrt  - f f-' L L ' , lU  S = X 

n N j r I.  c j . L ^ 1 
/ / *)  r c '■  X f,  / A 

/ / '■  V S I K u ■> 


I. = T'  X 

//S  T L.^-z  £ XtC  PLL  I CL  -•  .FC‘U  = • H S T f MAP*  ,'XEG1CN  = 1 ?:k 

//M'^l.SYSl^  cr  * 


C 


1 

2 

3 


C 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r. 

c 

r. 

c 

c 

c 

c 

c 


-r>-«L:r  Y.ii,/,L : ^^lottck  crcc- 

i:  I N ; C S It,  \ .11  S , 5 L 0 ) , F ( ‘i , j ) , c ( 5 , b DC  ) , r ( 5 , 5 c I 
!<  E w 1 1\  i:  i ' 
i:c  3 ! = 1 , 1 

tc  : 

f i £ A L ( I ■ , l ) 0(1,  n , L ( 1 , J ) , c ( I , .) ) , 0 ( 1 , J ) 

(-  r t<  N U ( A L 1 ■> . X,  ) 

cr  M I \ij‘ 

(.L\  r icc'^ 

C£Li  cs  I /:  ( -5.. , 1 1 , 1 12 1 ) 

CALL  SPPLLKA,'  •,  I , 7 . 5 , 5 . 0 , C , 5 , 5 2 T , I I 2 1 ) 

i.AlL  PLuT  ( ' ,-2.5,-  3 ) 

CALL  SPPLCMA,'  Sl.C,'  ',1,'  ' , 1 , 7.  5 ,5.0 ,0,  b ,bC'‘ , 1121  ) 

call  plC  M •■  . , -2  .S,  -3  ) 

CALL  5 P P I C T ( A , • • , I , f) , ' • , 1 , • ■ , 1 , 7 . 5 , 5 . : , C , b , b : P , 1 1 2 1 ) 

CALL  PLiiT  ( ..:,^S9) 

2 (LP 
t-fvC 

SLLXCUT  I Nl  SP(M  r.T  ( X , Xht'P  , K X , Y , YhDP  , K V , f HPP  , K F , 

1 A L fc  C X . A L e Y , L 1 A T YP  , N RC  W , NCCL  , JP£  N) 


SPPLCI  IC  A SUlPCLTINL  ThAf  WILL  PLCT  N'RO,s  NLN'PER  CF 
LI\Fj  PfP  Gw'iPI  . (V\R1TT£;-J  L'Y  GLP.PY  A.  LINDHCl.P) 

IFF  APGcPFMC  CF  TFL-  LI  21  Aft  FI: 

' IAm:  LC  TF  : ARRAY  CC  \ f A [ 1 WG  TFE  AESCISSA  TR  X 

X A L L f;  S . 

XH;R  A,\C  KX  CESGftlfE  the  CESlftEC  LEGEKC  FC^  TH^ 
x-AXf  f . 

xFCft  IS  n--  LITCPAL  EXPRESSIUN  CF  THE  LEGr\C  rc  3E 
px  i,%rr:  . 

iXX  IS  Ti-.  (\T£G£ft  CF  CHARACTERS  1C  XHCR  (LIMT 

CHAPACFiiRS  ) . 

Y HAHc  CF  THE  ARRAY  cCNTAIYirJG  ?h£  C.RCI\ATP  C«  Y 

Y Al  LE':  . 

VHCR  A.’JL  KY  CESCRMiE  THE  CCSIREC  LEGE'YC  FCR  THE 

Y - A X I S . 

YFCR  is  THE  literal  SXPRESSIGM  GF  THE  LEGEND  TC  CF 
PR  iCTi-L  . 

KY  IS  1HL  INIEGER  NUH‘*QR  GF  CHAftACURS  I\  YHCr  (LIMT 
2C  CF AftACrr <S  ) . 

FFCft  AM.  KF  C^SCRIhE-  TE'E  CtSlREC  LF.GE.NC  FCR  TF  C u«APM  HEADING, 

Li'x[T  iCv  characters. 

FFL'R  IS  (HE  literal  EXPRESSION  CF  fHF  HFAriNG  TC'  PE  PRINTED. 

KF  IS  TEE  INTEGER  NUMHER  GF  CHAKAIERS  IN  THE  HEADING. 

AI.CNX  I2  FFF  desired  LENGTH  CF  THE  X-AXiS  IN  INCHES  JLT  TL 
-XCEi-r.  (NGI'ES. 

'.DENY  IS  T,R  DESIREC  LENGTH  UF  TER  Y-AXIS  IN  INCHES  \CT 
M EXCEED  x.p  IncfES. 

LCIIYP  CCT^RVINES  THE  TYPE  CF  LINi.  YuG  GET. 

If  LI’  TYP  GREATER  THAN  C,  YCL  GET  A STRAIGHT  I PM 
CC  NNL-C  r I sG  A SYPfCL  PLT  CGKN  AT  E'/ERY  LINTYPTH 


r>  o o n o o r»  r>  o 
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I 


c :n  r.  I . 

ir  l:myp  wckf  cc^lal  if  2f  vr,i  wr:iLi;  i,l  t a 

Lir.'"'  CC  f I NC  LVriJY  SCCFNC  r-fllM. 

IF  LIMY.J  rtE^r:  FlJUAL  T L , YOU  .-.111,  GFT  ^ I.  irv<- 
(.M  Y . 

If  LI.MYP  LESS  IHA.'g  yGL  YsDIjLC  GET  A SyF'''U. 

LVf  RY  I IMYl'Tf-  PC  I\T  . 

.>pr.rt  IS  imemC-r  \u^cek  ihat  tells  the  rr 

li\es  i.rav.a  r:\  n-E  gra(>E'. 

'.CEL  IS  AK  INTEGtR  S'L^’LEP  M.P  Tf  t IvURhER  ff  GAIA  Rf  l jIS 
C\  EACL  LI\E. 

C JPE\  IS  Tf-u  TYPE  GF  P i m LSCIj  T (J  RakE  TEE  GRAPH, 

C I.".  PALL  PCIM  PhiJ, BLACK  R A P I f (,  C K A P fi  P = G,  i-TC.  ( R f -R 

C rC  TfiE  PLGITER  ^ANLAL  SLPRLIlI.  HY  TeE  CCPPLIFR  C P T E -?  LGCr-R 

C S I / E ) . 

C 

C 

LI^lNSIIA  XI\RC»v,\CLL)  , Y(NRCkN,:jC(jL  ) , XHDRI2:), 

1 YE  LR  ( ?C  ) , EE!;»  ( hA  ) 

CALL  S ERPT.L  I 2 . 5,  1 , . 1 5 , F HLK  , C . , K f ) 

call  PLCI  (2.  j,2.‘>  + ALr''fV,  ■)) 
gall  PLC  r ( 2 . 5 y ale  NX  ,2.  S+ AL  EivY  ,2  ) 

CALI  PLCT I 2 . 5 +ALENX ,2 .5 , r ) 
v,ALL  PLCT  I 2 . G , 2 .S  J ) 

call  SCALE2(X,ALE\X,rvRCU,NCCL,F'<X,CLX) 

CALL  SCALED  ( Y , Al_  e\Y  ,\RC»>  ,\CQL  ,fPY  , C LY  ) 

CALL  AXIS!..'  ,C.:,XECR,-KX,AL'-NX,.:.^,F-<X,LIX) 

L«LL  AXlC(^.:,L.:,YEL)R,KY,ALEGY,'TO.G,fPY,'  LV) 

(C  32  J=l,\Pf'K 

IFIJ  .\E.  1)  call  PLfJT(G.:,C.C,3) 

I N T • C = J 
I P C-  , = T 
ICCLC  = -1 
M = I APS ( L I NT YP  ) 

IF  ( ! I\T  YP  ) V , . 7 

t rvT^i 

7 r,r  = i 

K A = r 
K K = i 

IF  ( L I\ r Y P ) 11,12,13 

11  I P E A = ' 

ICCT.EA  = -i 
LSE»  = 1 

CC  TL  15 

12  AAr  ,CCL 

13  I PEG  A = 2 
ICCCEA  = -2 

lsy>  = : 

15  LL  3C  i = i,m:ci 

X\=(X|J,I)-F«X)/CLX 
Y\=( Y( J,  1 )-ERY  l/CLY 
IF  INA-M  ) 2;,  4 1,2  2 

?:  IF  (LSV»)  2 3,22,23 

21  call  SYPPCL  ( XN,YN,C.Cb, INTEG.C.O,  ICCLF) 

^A  - 1 

C C It,  2 5 

2 2 CALL  PI  r.  T 1 XN  , Y.\  , I PE'f!  ) 

2 3 NA  = NA  ♦ 1 

25  NF  = SFfKK 


4 


1 . 
J I 


/ 1 


it 


I P L = 1 r*  r 
I 1 M.L- 
CTM  J \Lt 


\u 


fN 


Pt  fL- 
l-NC 

SLHKf I r I \ 


/W  •J  \ Y 


a <L  i 


SCALC<>  ( 4PRAY  , 'i'CL  t^  ,^  <Cr.  , NCCL,  F I PST  ,Ci:'L  ) 

Cr  APPAY  CC.NTArNlRC  VALLLS  Tf  PC 
SCALF C . 

INOES  CvFP  XFloH  APPAY  IS  Tf; 


APT  ‘ 

1 v: 


1 I Y E ^ S I L \ 

S A V . ( i I = 1 . . 
SAv‘‘(?)  = ^.:, 

S A V . ( M = A . „ 

S A V L (A  ) = S . . 

s A ■-(-)  = p . : 

SAv:  (6  )=  1.  . 

S A vM  7 ) = .■’■:. 

F A r:  - c . :■ ! 

I N C ^ 1 
•\  p T S = ':c  Q L 
K = 1 *.i;S  ( I .(-  ) 

A =\'M  S'^N 

Y i = I . : f 6 

YF  = - 1 . 

CL  -'t  J J = 1 , 

Y ■-  =>>'<PA  Y ( J J 

Y \ = Y 

lL  t'f  I = 1 , - , K 
YS  = APAAy I o j , I 
IF  I Y : - Y S ) 

Y : = Y s 

cc  Tr  p-j 

IF  (YS-YM 

Y A = Y 

I F I Y V . : I 
1 F ( Y . . 1.  T 

CCAI  lALt 

Y : = Y 1 
YA  = YF 
F I R S 1 V = Y L 
IF  ( Yv,  ) 

F AL=:F  AL-1 
L EL  r Av= ( Y 


LCACrt-  I\ 

-Jii  SC  alec 

N'.Vt’CP  CF 

IaCRE^'EAT 
Pr lATS. 

tt  p < a Y ; A p r R , \ 


PLiNTS  r.j  p>t  sc^aEC. 

UF  LLCATIC'J  CF  SUCCESSIVE 


;0L)  ,SAVC-(  7) 


I ) 


) 

2? 


22,21 


2 S , 2 A 


Y f ) Y F = Y \ 

Y 1 ) Y i=y: 


I A , 1 i 


F IPSTV)/AXLF\ 


IF  (LELIAV)  7C,7C,AC 


a:  I =al:"g  1 ' ( cel  T AV  ) ♦ ICCC. : 


I 


[ FL  1 Av  = LEL  T Av/P-C 
LC  AS  I=1,L 
I S=  I 

IF  (SAVE!  I )-CEL  lAV)  AS,SC',5C 

ccAT  ikl: 

r E L i AV  = SAV:  (IS)  <’P 

I-  I rv=r ,:l  A i.s  T ( V.  /n'":L  r Av  + F AD ) 

1 zF  IPSTV  + ( uxL  :A+..  ..  D^DELTAV 
IF  (T-YSi)  S‘=,S7,S7 
F I P J T V^R-i  IJT  ( Y :/P  *FAC  ) 


J 


rv+(*vxl.  rNt.Cl  )<'CFLTAV 
If  ( r-y  I ) I ,^.7 

F 1 s = 1 S ♦ 1 
LC  rc  s: 

■j  7 f 1 S T V=r-  I f.  T V - A 1 \ T ( ( A <L  !.  \ + ( f 1 >7  S I V - V ,\  ) / ^ f L I A V ) / ? 
1 1:  E L I A V 

IF  ( Y r ( -i  i I i/  ) 5 f!  f •■)  r * “j  y 

-jh  f I4S  rv=G  . : 

5-y  IF  ( I\C  ) 7.  i , M ,f. 

O 1 f Ik':,  rV=F  IF  S r y t A [,NT  ( AXL  F:.'i  + . 5 ) ‘J'ClL  TA  y 

i:  E L r A V = - n F L T A V 

0.  r IkbT-F  I I y 
C E L = f;  F LI  A V 
6 7 k t T L k \ 

7C  LEL  r AV  = ^' .1.:  *F  I kb  T'y 

Lf:l  I Ay  =A'<S  ( L'EL  I A v/AXLK  G ) + L . 
cc  rc  AC 
1 1\  c 

//CC.PLCr,.  1 E[  Lb\AfV  = kECAF.(,AIEk, 

//  r I SP=  ( itV.  ,C  AT  Lb  ) ,LM  r^PLIHC  I SK  , 

//  LF  F--PL»<  6 I / £ = , SPACI:  = ( C Yl  , E , kf  St  ) r I A p I I = k f-  T P P = P 

//GL.F  Tl'Jf  Cwl  rc  CSN--FKCA  P.  P IKGGtL  I Sk  = EHK  ,LM  1=^T  API 

/ /c(,. SVC  I . nr  y 


APPENDIX  II 


VERIFICATION  OF  THE  COMPUTER  CODE 


The  following  test  cases  were  used  to  show  that  the  Thrust  Malalign- 
ment Computer  Code  will  give  accurate  results  in  the  prediction  of  motion 
about  a fixed  point.  Each  case  was  solved  by  an  exact  analytical  techni- 
que and  also  by  the  computer.  The  results  of  each  method  of  solution 
are  indicated. 


Test-1 


1 

i 


1 


Calculation  of  the  Period  of  Motion  of  a Simple  Pendulum  in 

a Gravitational  Field 

Consider  a simple  pendulum  as  shown  in  Figure  29.  Mass  m is  a fixed 
distance  1 from  the  pivot  point  0 in  coordinate  system  xyz.  XYZ  is  the 
inertia  reference  frame.  0^  and  ii;^  are  the  original  nutation  and  pre- 
cession angles.  The  pendulum  is  released  in  a gravitation  field  with 
an  acceleration  vector  in  the  -Z  direction.  1 = The  problem 

is  to  study  e[T]  using  the  computer  program  for  a single  case  of  0^, 

The  period  of  motion  for  the  pendulum  will  be  compared  to  the 
standard  solution  for  the  problem. 

Standard  Solution 

The  exact  solution  may  be  obtained  as  follows  [3]: 


complete  elliptic  integral  of  the 
first  kind. 


k^-  sin^  <}i 


T = 4 


Vt 


K(k) 


where. 


K(k)  = F j J,  kj 


k = sin 


K(k) 


/ 


T 

tt/2 


V 


1 ■ 


t 60 

L. 


Figure  29.  A Simple  Pendulum  Os 


r 
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Consider  a pendulum  with  the  following  characteristics: 
g = 32.1  ft/sec2 
m = 1.0  slug 

w = 32.1  lbs  I = .0001  slug-ft2 

1 = 1 .0  ft  r = 1 .0  slug-ft2 


^0  = 


0 = 80.0' 
0 


K(k)  = F (|j,  sin  40°  = 1.7868 


= 

^ V 32.1 


1 .7868  = 1 .2615  sec. 


The  final  results  are: 


1.2618  sec 

computer 


W 


Test-2 


i 


[ 

1 


I 


Spinning  Gyroscope  in  a Gravitational  Field 

Consider  the  problem  of  determining  the  motion  of  a spinning  gyro- 
scope in  a gravitational  field  as  shown  in  Figure  30.  The  gyroscope  is 
a 1 in.  thick,  6 in.  diameter  steel  disk  mounted  midway  along  a 12  in. 
spindle.  The  spindle's  spin  axis  is  always  coincident  with  the  z axis 
in  xyz.  The  inertia  reference  is  the  XYZ  coordinate  set. 

Find: 

the  rate  of  precession, 
the  nutation  amplitude, 
the  nutation  period,  T. 

where, 

6q  = 7500  rpm  = 785.398  rad/sec 

.1^  = 0° 

^0 

e = 15° 

0 

1 = .5  ft 
r = .25  ft 
h = .0833  ft. 
m = 2 slugs 
g = 32.16  ft/sec2 

I = ^=  .0625  slug-ft2 

r = nL_(^^  = .5324  siug-n2 

w = 64.32  lbs 
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The  results  of  the  computer  solution  compare  with  the  exact  soluti 
as  follows: 


Nutation  Period 
Nutation  Amplitude 
Precession  Rate 

*SEE  REFERENCE  [4]. 


COMPUTED  SOLUTION 
.069  sec/cycle 
.00378  radians 
.642  rad/ sec 


EXACT  SOLUTION* 
.068  sec/cycle 
.00374  radians 
.655  rad/sec 


REFERENCES 


1.  Shames,  I.  H.,  Engineering  Mechanics -Dynamics.  Prentice-Hall,  Inc., 
Englewood  Cliffs,  New  Jers^,  (1966).~ 

2.  Meriam,  J.  L.,  Dynamics-2nd  Edition.  John  Wiley  and  Sons,  Inc., 

New  York,  (1971] . 

3.  Greenwood,  D.  T.,  Principles  of  Dynamics,  Prentice-Hall,  Inc., 
Englewood  Cliffs,  New  Jersey,  1[^T965) . 

4.  Scarborough,  J.  B.,  The  Gyroscope,  Theory  and  Applications, 

Interscience  Publishers,  New  York,  (1968).  ~ 


66 


